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§ : ABSTRACT 

We investigate a corona coupled with underlying disk through magnetic 
field and radiation field, and present emergent spectra. Due to buoyancy the 
magnetic flux loop emerges from the disk and reconnects with other loops 
in the corona, thereby releasing the magnetic energy to heat the coronal 
(N| \ plasma. The energy is then radiated away through Compton scattering. 

By studying the energy balance in the corona, transition layer and disk, we 
determine the fraction (/) of accretion energy dissipated in the corona for 
^ \ given black-hole mass and accretion rate, and then determine the coronal 

and disk variables. This allows us to calculate emergent spectra through 
Monte Carlo simulations. The spectra are then determined as functions 
of black-hole mass and accretion rate. We find two types of solutions cor- 
responding for hard spectrum and soft spectrum. In the hard-spectrum 
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solution, the accretion energy is dominantly dissipated in the corona, sup- 
porting a strong corona above a cool disk; The hard X-ray spectral indices 
are the same for different accretion rate, i.e. a ~ 1.1 (F u oc v~ a ). In 
the soft-spectrum solution, the accretion energy is mainly dissipated in the 
disk. The coronal temperature and density are quite low. Consequently, the 
spectra are dominated by the disk radiation peaked in UV and soft X-rays. 
For low-luminosity systems (L < 0.2LEdd), there exists only the solution of 
hard spectra; While for high-luminosity systems (L > 0.8ZvEdd), there ex- 
ist both solutions of hard and soft spectra. For middle-luminosity systems 
(0.2L Edd < L < 0.8L Edd ), besides the hard spectra, moderately soft spectra 
composed of an inner soft-spectrum solution and an outer hard-spectrum 
solution may occur, the softness of which increases with increasing luminos- 
ity. The hard spectra are close to the observed spectra in Seyfert galaxies 
and radio-quiet QSOs. The composite spectra may account for the diversity 
of broad band spectra observed in narrow-line Seyfert 1 galaxies. 



Subject headings: accretion, accretion disks-galaxies: nuclei-X-rays:galaxies 
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1. Introduction 

It is commonly believed that the hard X-ray radiation of active galactic nuclei 
(AGNs) or galactic black hole candidates (GBHCs) arises from hot gas around accret- 
ing black holes, either in forms of hot accretion flow or accretion- disk corona (e.g. Liang 
& Nolan 1984; Mushotzsky, Done, & Pounds 1993; Narayan, Mahadevan, & Quataert 
1998). The structure of hot accretion flows, such as advection-dominated accretion flow 
(ADAF), has been extensively investigated based on the vertically one- zone approxi- 
mation (Kato, Fukue, & Mineshige 1998 and references therein). The ADAF model 
(Ichimaru 1977; Narayan, Yi & Mahadevan 1995; see also Blandford & Begelman 1999; 
Igumenshchev & Abramowicz 2000 for modifications) can reasonably fit the observed 
spectra of low luminosity AGNs and X-ray novae in quiescence. Nevertheless, spectral 
features observed in ordinary AGNs, the big blue bump, the soft X-ray excess and hard 
X-ray continuum, and the 6.4 keV fluorescent iron lines, show strong observational ev- 
idence for hot gas coexisting with cool gas in the vicinity of accreting black hole (e.g., 
Mushotzky et al. 1993), which is presumably a corona lying above an accretion disk. 
In disk-corona models, the cold disk is embedded in the hot corona in a plane-parallel 
slab. Except for the ion-illuminated disk (Deufel & Spruit 2000; Deufel, Dullemond, & 
Spruit 2002), previous disk corona models for AGNs or GBHCs (Haardt & Maraschi 
1991; 1993; Nakamura & Osaki 1993; Svensson & Zdziarski 1994; Dove et al. 1997; 
Kawaguchi, Shimura, & Mineshige 2001) need to assume a large fraction of accretion 
energy to be released in the corona though the details of the coronal heating mechanism 
remain unclear. 

Detailed investigation on the interaction between an accretion disk and a friction- 
heated corona (Liu et al. 2002b) shows that the energy deposit in the corona is not 
enough to keep itself hot above the disk against strong Compton cooling in AGN 
systems. Other heating besides the viscous one is required to maintain such a hot 
corona lying above the disk and to produce the observed X-ray luminosity. Comparison 
of the thermal energy capacity in the corona and observed power of AGNs also shows 
energy shortage in the corona (Merloni & Fabian 2001). One promising mechanism to 
heat the corona may be the magnetic reconnection (e.g. Di Matteo 1998; Di Matteo, 
Celotti, & Fabian 1999; Miller & Stone 2000; Machida, Hayashi k Matsumoto 2000) as 
a result of magnetorotational instabilities in the disk and the buoyancy of the magnetic 
field (Tout & Pringle 1992; Miller & Stone 2000). Despite of the complex radiation 
and energy interaction between the disk main body and the corona, magnetic fields 
seem to play essential role in producing time variations and spatial inhomogeneity (e.g. 
Kawaguchi et al. 2000). 

In a recent study (Liu, Mineshige, & Shibata 2002a, hereafter Paper I), we show 
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that, the magnetic reconnection can indeed heat the corona to a temperature around 
10 9 K and produce the observed X-ray luminosity. In the present paper, we describe 
the model in a more consistent way (Sect. 2), and then calculate the spectra from 
the corona and disk by Monte Carlo simulations. The computational results on the 
coronal properties and emergent spectra are presented in Sect. 3. Discussion on spectra 
is presented in Sect. 4 and our conclusions in Sect. 5. 



2. The model 

In our disk-corona model (Paper I), the disk is assumed to be a classical Shakura 
& Sunyaev (1973) disk. The corona is a plane-parallel corona which is tightly coupled 
with the underlying disk. In the disk, magnetic fields are continuously generated by 
dynamo action. Owing to the buoyancy magnetic flux loops emerge into the corona 
and reconnect with other loops. Thereby, the magnetic energy carried from the disk is 
released in the corona as thermal energy. If the density of corona is not high enough, 
heat is conducted by electrons from the corona to the chromosphere, resulting in mass 
evaporation. Once the density of corona reaches a certain value, Compton scattering 
becomes dominant in cooling, and eventually an equilibrium is established between 
the reconnection heating and Compton cooling, a stationary corona is built up. The 
equations describing these processes in the corona and at the interface of disk and 
corona are (see Paper I), 
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In stead of r in Paper I, we introduce effective optical depth, r*, which includes the 
isotropic incident photons undergoing up-scattering in a plane-parallel corona, 

t* = \ T r = \ T nax(- (3) 

with A T being order of unit. 

Eq.(l) and Eq.(2) determine the temperature and density in the corona as function 
of energy densities of magnetic field and radiation field. 

The magnetic field strength B is derived by assuming equipartition of gas energy 
and magnetic energy in the disk, 

(3 = ndisk B fc 2 Tdisk ~ 1. (4) 
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For a radiation pressure-dominated disk, 

B = 1.47 x 10 3 a 1 [m o .i0 (1 - r&G; (5) 

For a gas pressure-dominated disk, 

B = 7.18 x 10 4 a~f°(3; h m^ [m ai 0(l - f)]K^G, (6) 

where / is the energy fraction dissipated in the corona; = 1 — \J~^ and is taken 
to be the last stable orbit 3i?s; m 8 , ^o.ij r i0) tt o.i> A are the mass of black hole, the 
accretion rate, the distance, the viscous coefficient, and the equipartition factor in units 
of 1O 8 M , 0.1M Edd , lOi^s, 0.1, and 1, respectively. 

The energy density of soft photon field, C/ ra( j, is contributed from both intrinsic 
disk radiation U^ d and reprocessed radiation U™d °f backward Compton emission from 
the corona 

in 4 4 3GMM(l-/)0 

^rad - aJ eff ~ C " g rf ^ ^ ( 7 ) 

= 1.14 x 10 5 m 8 1 m o .i0(l — /)r 10 3 ergs cm 3 , 



UH d = 0A\ u Ub- (8) 

In Eq.(8), a new factor X u is added to the evaluation of seed field in Haardt & Maraschi 
(1991; 1993), which includes two effects on the soft photon energy, i.e. the dependence 
on the coronal temperature and optical depth, and the ratio of Alfven speed to the 
light speed since the magnetic reconnection releases energy at Alfven speed while this 
energy is radiated away at light speed. 

Now we are allowed to derive two sets of coronal solutions, corresponding to a 
radiation pressure-dominated disk and a gas pressure-dominated disk, respectively. 
From Eqs.(l), (2), (5), and (7), we derive a coronal solution above a radiation pressure- 
dominated disk, 

Tx = 1.36 x lOVfA ~ 8 Karris « 2 [m ai 0(l - f)]~ r?^l K, (9) 

_1B _3 _1 _ 19 75 _3 

rai = 2.01 x 10 9 a rPi 3 A r 5 m 8 16 [m 0A (f)(l - f)}' 2 rf Vcm" 3 , (10) 

where the length of the magnetic loop £ 10 is in unit of 10i?s- 

From Eqs.(l), (2), (6), and (8), we get a coronal solution above a gas pressure- 
dominated disk, 

T 2 = 4.86 x 10 9 a~ff3;h;hu l4 mf [m o .i0(l - f)]™ r^l\ Q K, (11) 
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n 2 = 2.55 x lO 10 ^ 40 /?! ~ 4 \ T '\u~ 2 m s 40 [m o .i0(l - /)]* V° Vcm (12) 

The fraction of gravitational energy dissipated in the corona can be deduced from 
its definition (Paper I), 
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Replacing B and n by Eqs.(5) and (10), we obtain an equation concerning / for a 
radiation pressure-dominated disk, 

_45 _9 1 __9_ 225 3 

/ = 1.45(1 - /)- 2 a .f ^ «A 4 m 8 32 (mo.i0)-Vo 4 4) = Ci(l - /)~ 2 . (14) 

Replacing 5 and n by Eqs.(6) and (12), we get an equation for a gas pressure-dominated 
disk, 

11 _93 -I! I I I! 1 _JLL 3 

/ = 4.70 x 10 4 (1 - /) "c^fft 8 A 4 A 4 m 8 80 (m o .i0)^r lo 16O £f o . (15) 

Eq.(15) has solutions for any accretion rate; while Eq.(14) has solutions only if c\ < 
4/27, which gives, 

M — M _3 J_ _JL 15 1 

-r— > 0.21a ? *\?m 8 32 fVS4. (16) 

This is to say, at high accretion rate, the corona can coexist with a radiation pressure- 
dominated disk or a gas pressure-dominated disk; while at low accretion rate, the 
corona can coexist only with the gas pressure-dominated disk. The value of critical 
accretion rate is related to the distances. For a = 0.3, (5 — 1, A T = 2.7 (see next 
section), and M = 1O 8 M , a radiation pressure-dominated disk exists around distance 
R ~ QR S if M w 0.3M Edd ; extends to 50i2 s if M w 1.2M Edd . We will see in the 
following section that the spectra from a radiation pressure-dominated disk + corona 
and from a gas pressure-dominated disk + corona are of much difference. 

Solving Eq.(14) (or Eq.(15)) for / for a given black hole mass and accretion rate, 
we determine the coronal quantities from Eqs.(9) and (10) (or Eqs.(ll) and (12)). The 
disk quantities can also be determined. We are now ready to calculate the spectrum. 



3. Spectra calculated from Monte Carlo simulations 
3.1. The method of Monte Carlo simulations 



We use the Monte Carlo simulations to calculate the emergent spectrum. Our 
model consists of a cold disk, which produces blackbody radiation at each radius with 
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temperature Tr, 

u; \! ;/ I - f) r ( , 



4 3GMMt(l-f) c / 3GMM0(1-/) c , 



The hot corona lying above the disk is approximated as an optically thin slab with defi- 
nite electron temperature and scattering optical depth for each radius. The blackbody 
photons emitted from the cold disk go into the corona, some of them pass through 
without being scattered, some of them are Compton up-scattered. By tracing the pho- 
ton's motion in the corona, we separately record the photons emerging from upper and 
lower boundaries of the slab corona: the former are the ones to be observed and the 
latter impinge the disk and are reprocessed as blackbody emissions. Since the accretion 
disk is sufficiently dense, the albedo is neglected in our simulations. The spectra from 
different radii are added up as the total emergent spectrum. 

The method of the Monte Carlo simulation is based on Pozdnyakov, Sobol' & 
Sunyaev (1978). We restrict our consideration to a thermal corona where the elec- 
trons have a Maxwellian distribution, and the Compton scattering is the main cooling 
process. Since the electron-scattering optical depth in the corona is less than one, we 
introduce the weight w as described by Pozdnyakov et al. (1978) in order to efficiently 
calculate the effects of multiple scattering. We first set wq — 1 for a given soft photon, 
then calculate the escape probability P of passing through the slab. The quantity of 
WqPq are the transmitted portion and is recorded to calculate the penetrated spectrum 
or reprocessed photons according to the escape direction of the photon. The remaining 
weight w\ = wo(l — Pq) is the portion that undergoes at least one scattering. If we 
write the escape probability after the n-th scattering as P n , the quantity w n P n is the 
transmitted portion of photons after the n-th scattering, and is recorded as upward 
or downward transmitted spectrum. The remaining portion w n (l — P n ) undergoes 
the (n + l)-th scattering. This calculation is continued until the weight w becomes 
sufficiently small. The whole process is simulated by the Monte Carlo method. 



3.2. Iterative computations for consistency 

Once the coronal temperature, the optical depth and the seed soft photons are 
given, the Mento Carlo simulation can be performed and the Compton scattering spec- 
trum and luminosity coming out from both upside and downside of the corona are 
obtained. The questions arise, does the upward-outcoming luminosity equal to the 
gravitation energy released by the accretion? does the energy density of downward- 
outgoing photon plus intrinsic disk radiation equal to the presumed energy density 
of seed soft photons? If the seed photon energy and the effective optical depth are 
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correctly chosen in deriving the corona temperature, density and optical depth, the 
answers should be affirmative. A consistent model should correctly describe the nature 
of corona, i.e., the temperature, density, optical depth, etc. so that the simulations 
give consistent results. 

In Paper I, we did not introduce A r and X u , which implies X T — 1 and \ u = 
1. This is to say, it is assumed that the reprocessed seed field has energy density 
^rad = 0.45 2 /87r and the effective optical depth is the vertical scattering optical depth 
t = tiutL However, unlike in a spheric geometry, in a plane geometry the actual 
optical depth (we called it the effective optical depth) should be larger than the ver- 
tical one, since the incident photons are isotropic and undergo longer path than I if 
their directions are not normal to the corona plane. Furthermore, A r is expected to be 
independent on the system parameters (e.g. accretion rate) since it is only a geomet- 
ric effect. The energy density of seed photon field contributed from the reprocessed 
photons, U™ d = 0AB 2 /8n, is estimated for specific range of parameters (e.g. Haardt 
& Maraschi 1993). The precise value depends on the temperature and optical depth. 
In order to get consistent results, we need to adjust the presumed energy density of 
soft photons and the effective optical depth of the scattering medium. This is why we 
introduce the coefficients A r and \ u in this study. For simplicity, we set A M the same 
value at all distances though the temperature and density of the corona change a little 
along the distance as shown in Paper I. This is reasonable since the reprocessed photons 
are not required to get out of the disk at the same distance where they impinge. 

For a given mass of black hole and accretion rate, we first perform the computations 
for the case of gas pressure-dominated disk. As the first step, we start the computation 
with initial values A r = 2 and X u — 1. For one grid point of radius R, we solve Eq.(15) 
numerically and obtain a solution for /; then calculate the corresponding temperature 
T and density in the corona from Eqs.(ll) and (12) and hence the scattering optical 
depth r, and the radiation temperature Tr of seed photon field by combining Eqs.(8) 
and (6). With T, r, and Tr, we perform the Monte Carlo simulation, record the local 
upwards luminosity and downwards luminosity. Then we change R to the next grid 
point and repeat the calculation for /, then for T, r and Tr, and then Monte Carlo 
simulation. Finally all the local outgoing photons are added up and the integrated 
upward luminosity L up and downward luminosity Ldown are obtained. Meanwhile, the 
integrated soft photon energy L soft = J^ s 2nRaT^dR and the released gravitational 
energy L G = J^ s 2irR(3GMM(f)/8irR 3 )dR are also calculated. The second step is to 
check the consistency, i.e. Ldown ~ -^soft an d L up LqI If yes, we find the correct A T 
and \ u and hence the consistent solution for the corona, and stop the computation; 
if no, set new values for the (n + l)-th A T and \ u from the n-th results: A Ui „ + i = 
(-^down,n/^soft,n)A„,„ and A T>n+ i = (L up>n /L G )\ T>n , then repeat the computation from 
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the first step until the consistent conditions L doW n R3 L so f t and L up m Lq are fulfilled. 

In the case of radiation pressure-dominated disk, We solve / from Eq.(14) and 
then the temperature and density from Eqs.(9) and Eqs.(lO). The iteration procedure 
is similar to the case of gas pressure-dominated disk. The difference is in that we 
only need to check whether L up « L G is fulfilled or not, since the contribution to the 
seed photons from the reprocessed photons are not so important as the intrinsic disk 
radiations. Thus, \ u is simply set to unity and is not adjusted during the iteration. 

3.3. The self-consistent coronal structure 

For given parameters, M = 1O 8 M , a = 0.3, £ = 10i?s, P — 1, we perform 
computations for a series of accretion rates. As we discussed in Sect. 2, computations 
show that, for low accretion rates (M < 0.3M Edd ) there is only gas pressure-dominated 
solution; for high mass accretion rates there are both gas pressure-dominated and 
radiation pressure-dominated solutions (say, M = 2MEdd); for middle accretion rates 
(say, M = 0.5MEdd), the gas pressure-dominated solution exists at all distances, while 
the radiation pressure-dominated solution only exists at small distances. As an example 
for low accretion rate system, Fig.l shows how much accretion energy is dissipated in 
the corona and in the disk for M = 0.1M E dd- It is shown that almost all the accretion 
energy is dissipated in the corona; less than one thousandth is dissipated in the disk. 
The corresponding coronal structure is plotted in Fig.3. The temperature is ~ 10 9 K 
and the density ~ 10 9 cm -3 , which are a little lower than that in Paper I since here 
we take into account the effective optical depth and consistent seed photon field. The 
effect optical depth is around 0.7 and Compton y-parameter 0.6. 

As an example for high accretion rate, Fig. 2 and Fig. 4 show the energy fraction 
and coronal structure for an accretion rate of 2Mgdd- For a gas pressure-dominated 
disk, the coronal temperature and density are almost the same as that for low accretion 
rates. Only the magnetic field and Alfven speed increase with accretion rates. This 
feature is the nature of gas-pressure solutions because / depends on M very weakly 
(see Eq.(15)), so are the temperature (see Eq.(ll)) and density (see Eq.(12)). In other 
words, the coronal structures with a gas pressure-dominated disk underneath are quite 
uniform. Higher accretion rate only results in faster energy transfer and dissipation in 
the corona. Such features can be understood as follows. A corona coupled with the 
disk through the magnetic field and radiation field can self-adjust to an equilibrium at 
a certain temperature. The magnetic field transports almost all the accretion energy 
to the corona and releases the energy there through reconnection. At the beginning 
of reconnection, the electrons are suddenly heated to a very high temperature. Even 
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though there is no many seed photons to be scattered, the Compton cooling power is 
still quite strong due to the very high temperature. Thus, the seed photon's energy is 
greatly amplified. Some fraction of the scattered photons get out of the corona from 
the backside and are reprocessed in the disk as blackbody radiation. Then the seed 
photon field is strengthened and the electron temperature decreases due to efficient 
Compton cooling. The equilibrium is finally reached when the heating by the mag- 
netic reconnection balances the Compton cooling. Such an approach to the equilibrium 
can be clearly seen from our iterative computations: when \ u is set to be a small value, 
the temperature derived from the model is high. With this high T we then obtain high 
backward-scattered photon luminosity (by the Monte Carlo simulation). Taking this 
high luminosity as the seed soft photon luminosity, the newly derived T becomes low 
and hence the soft-photon field contributed from the backward Compton radiation is 
closer to the true value. By repeating this procedure, the backward-scattered luminos- 
ity quickly converges to the seed soft luminosity, and a stationary state is reached. 

In table 1 we list for a gas pressure-dominated disk the coefficients A u , A T , the 
emergent luminosity L up , the seed soft photon luminosity L S oft, and the amplification 
factor A = (L up + L down ) / L soft , where L down is the downward/backward Compton 
luminosity. We find that the Compton amplification factors for different accretion 
rate are almost the same, A ~ 2.3. Soft photons are scattered in the corona, gain 
energy by a factor of A — 1 ~ 1.3 and escape from both sides of the corona. The 
upward-escaped photons carry a little more energy (L up ) than the downward-escaped 
photons (ps Lsoft), the latter is reprocessed as seed photons. The coefficient A M for the 
seed photon field increases with increasing accretion rate. Comparing X u and Alfven 
speed, we find that \ u /Va does not change with M. This indicates that the reprocessed 
photon field is actually proportional to the Ub(Va/c), the factor Va/c is due to the lag 
of magnetic heating (at a speed of Va) to radiation cooling (at a speed of c). We also 
find A r ~ 2.7, independent on M. This is to say, the effective optical depth caused by 
the slab geometry is as 2.7 times large as the vertical optical depth. 

In contrast, for a radiation pressure-dominated disk, a large fraction of the accre- 
tion energy dissipates in the disk (see Fig. 2), the corona is quite weak with temperature 
being a few 10 8 -K" and density ~ 10 8 cm~ 3 (see Fig. 4). The seed photons are mainly 
from the intrinsic disk radiation. The disk and corona are not coupled at the same way 
as that for a gas pressure-dominated disk. In such a strong disk + weak corona system, 
the Compton scattering in the corona is quite weak, A T ps 1 and A m 1. We note that, 
the energy fraction dissipated in the corona decreases with increasing accretion rate, 
the temperature and density in the corona also decrease with M. This implies that the 
corona above a radiation pressure-dominated disk is weak at luminous systems. 
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Table 1: Computational results for a corona above a gas pressure-dominated disk at 
different accretion rates (where MEdd = -^Edd/O.lc 2 ) 



0.01 2.35 1.10 2.74 4.43 x 10 43 3.30 x 10 43 

0.05 2.34 1.89 2.72 2.21 x 10 44 1.65 x 10 44 

0.1 2.34 2.38 2.71 4.43 x 10 44 3.30 x 10 44 

0.5 2.33 4.08 2.68 2.22 x 10 45 1.66 x 10 45 

1.0 2.34 5.12 2.68 4.45 x 10 45 3.30 x 10 45 

2.0 2.35 6.40 2.68 8.84 x 10 45 6.57 x 10 45 



3.4. The emergent spectrum 

In last section we show that there are two types of solution for the disk and corona 
structure. When the accretion disk is presumed to be gas pressure-dominated, we get 
a solution / ~ 1, indicating that most of the gravitational energy is transferred to 
the corona by strong magnetic field and fast reconnection. The disk is quite cool, 
with midplane temperature of a few 10 4 K. The corona is heated up to a temperature 
T ~ 10 9 K by the magnetic reconnection, and the density becomes 10 9 cm~ 3 through 
efficient mass evaporation at the interface between the disk and the corona. The 
dissipated energy in the corona is emitted away by Compton scattering. Part of the seed 
photons are upward scattered as the emergent spectrum, and part of them are scattered 
backward, which are reprocessed in the disk surface layer and are emitted as blackbody 
spectra. Here we assume that the irradiation by the backward Compton radiation does 
not largely influence the internal structure of a radiative disk but only changes the 
temperature of the surface layer (for details see Tuchman, Mineshige, & Wheeler 1990). 
Therefore, the magnetic energy determined by the equipartition in the disk midplane 
is not affected by the irradiation. Such a highly heated corona produces strong hard X- 
ray emission. The corona and disk are similar to the "two phases" discussed by Haardt 
and Maraschi (1993). We call this solution as hard state since its spectrum is hard. 
Fig. 5 shows the spectral energy distribution of hard state for different accretion rates 
corresponding to different luminosities. [Here we need to point out that the energy 
conversion coefficient rj by accretion is not exactly 0.1 in our case. For given accretion 
rate we integrate the gravitational energy by Lq = f™ Hs 2irR(3GMM(l)/87rR 3 )dR. The 
total luminosity from both side of the disk and corona should be 2L G . For R* = 3i?s, 
i] = 2L G /Mc 2 ~ 0.07. Thus, even for M = M Edd = L E dd/0.1c 2 , the luminosity is 
less than the Eddington luminosity] From the figure we see that for different M the 
spectral shapes in the hard state are very similar. The hard X-ray spectral index is 
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around 1.1, which is a little steeper than that observed in Seyfert galaxies. 

The "constant" spectral index is caused by the "uniform" corona. As we show 
in Sect. 3. 3, the coronal structure at hard state is nearly independent on the accretion 
rate and the Compton y-parameter (= (4/cT/m e c 2 )r) hardly changes with M. Our 
computational data show that the radial distributions of ^-parameter are indeed the 
same for different M; y reaches the maximum of y ~ 0.6 in the inner region and 
decreases with distances (see Fig.3). Therefore, the hard-state spectra have similar 
spectral shapes. 

On the other hand, if the disk is presumed to be radiation pressure-dominated, 
there exists no solution for / at low accretion rate. At high accretion rate, the solution 
for / is far less than 1, indicating that magnetic buoyancy does not transfer much energy 
from the disk to the corona, and there is only a weak corona above the disk with low 
temperature and low density compared to those in the hard state. The gravitational 
energy is dominantly released in the disk as multi-color blackbody radiation. These 
photons go through the very weak corona, almost all of them penetrate the corona 
without being scattered. The spectral energy {vL v ) peaks at around UV to soft X-rays, 
dropping steeply at high frequency. We call this solution as the soft state. Fig. 6 shows 
typical spectrum of the soft state for M = 2MEdd (corresponding to L — 1.4LEdd)- The 
soft-state spectrum is not much different from pure-disk spectrum. For comparison the 
hard-state spectrum for the same M is also shown in the figure. Obviously, the spectral 
distributions at soft and hard state are very much different. One may question whether 
the Bremsstrahlung emissions contribute more to the hard X-rays in this state. Our 
estimation shows that the Bremsstrahlung emission power is about 5 order of magnitude 
less than Compton radiation because the density in the soft-state corona is so low. 

Furthermore, there are composite solutions, in which the disk is composed of 
an inner radiation pressure-dominated region and an outer gas pressure-dominated 
region with a weak inner corona and a strong outer corona above. Whether such a 
discontinuous disk + corona can dynamically exist is to be studied. Nevertheless, 
the spectra from such a composite corona + disk show moderate UV-X-ray spectral 
indices. In Fig.7 we plot such a spectrum for M = 0.5MEdd (corresponding to L — 
0.35LEdd)- It is shown that the inner region produces soft spectra dominated by the 
disk radiation and the outer region produces typical hard-state spectra dominated 
by the corona radiation. Thus, the composed spectrum is moderate: its hard X-ray 
spectrum has the same spectral index as that of hard state; while the ratio of UV 
and X-ray luminosity, or rather, the UV-X-ray spectral index, is larger than that of 
hard state. With M increases, the inner radiation pressure-dominated region extends 
outwards and the outer gas-pressure dominated region shrinks. The overall spectrum 
then becomes softer. 
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As discussed in Sect. 3. 3, the hard-state solution can exist in both low- and high- 
luminosity objects; while a soft-state corona extending to 50-Rg can only appear at 
high accretion rate, M > 1.2M Edd (or L > 0.8L Edd ). Systems accreting at a rate 
0.3ME dd < M < 1.2MEdd ( or 0.2LE dd < L < 0.8Z-E dd ) can be in a state between hard 
and soft with UV-X-ray indices varying with accretion rates. Therefore, our model 
predicts that, two spectral states are possible for accretion rate above 1.2ME dd ; below 
0.3M E dd only hard state. These features of hard and soft states are summarized in 
Table 2. 



Table 2: The features of disk and corona in the hard and soft spectral states 



Spectral 


Range 


/ Pressure Typical coronal values 




a x 


state 


of L 


in disk T(K) n(cnr 3 ) r* 






Hard 


Always 


~ 1 P| isk ~ 10 9 ~ 10 9 0.6 


~ 1.4 


~ 1.1 


Soft 


> 0.8L Edd 


< 0.1 P disk Low, decrease with L 


> 2 




Moderate 


0.2-0.8L Edd 


outer region hard, inner region soft 


> 1.4 


~ 1.1 



The parameters in our model are the mass of black hole M, the accretion rate 
M, the thickness of corona £ (equivalent of averaged length of magnetic loops), and 
the equipartition coefficient (5. Above study shows the corona structure and spectra in 
dependence on M for given parameters M = 10 8 M Q and £ = 10i? s in the region from 
3i? s to 50i? s . From the equation concerning / and the solutions of T and r(oc nl), we 
know /, T and r and hence the spectra only very weakly depend on £, especially at 
the hard state. In Fig. 8 we plot the hard-state spectra for £ = lORs and £ = 20R S . 
The two curves overlap. Similarly, the spectra have weak dependence on the mass 
of black hole. But, considering large change of black- hole mass over a few orders of 
magnitude, the spectra may be somewhat different. In Fig. 8 an example of hard state 
for M = 10 5 M Q is shown. Compared to the case for M = 1O 8 M , the disk component 
shifts towards high frequency, whereas no obvious difference in the shape of hard X-ray 
spectra contributed by the corona. For a 10M Q black hole we expect similar hard 
X-ray spectra except for the low-frequency component which peaks at soft X-rays due 
to high effective temperature in the disk. The dependence on (3 seems large because 
it directly determines how much of the accretion energy transferred to the corona 
through magnetic field. Owing to poor knowledge on the equipartition coefficient, we 
don't discuss this issue here. 
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4. Discussion on the spectra 

The observed spectrum in X-ray 2-20 keV is close to a power law, with an averaged 
spectral index m 0.7 ± 0.15 for Seyfert galaxies (Mushotzky 1984; Turner & Pounds 
1989). Taking into account the deconvolution of the Compton reflection hump around 
^30keV (George & Fabian 1991; Williams et al. 1992; Mushotzski et al. 1993; Petrocci 
et al. 2000), the spectral index of the underlying power law for Seyfert galaxies is 0.9 
to 1.0 (e.g. Nandra & Pounds 1994; Pounds et al. 1990). The X-ray spectral index for 
radio-quiet QSOs is slightly higher, ax ~ 1-0 (e.g. George et al. 2000). The observed 
broad band spectral index measured between 2500A and 2keV is 1.25 for Seyferts 
(Walter & Fink 1993) and 1.5 for QSOs (Yuan et al. 1998). Our model shows that at 
hard state the X-ray spectral indices between 2-20keV are around 1.1. The spectral 
indices between the hump (around UV) and 2keV, a px , are around 1.4. The hard-state 
spectra are roughly in agreement with the observed spectra of Seyfert galaxies and 
QSOs. 

The composite spectrum, which occurs only at relative high accretion rates and 
becomes more and more close to the soft spectrum with increase of M, may explain 
the diversity of spectra observed in narrow-line Seyfert 1 galaxies (NLSls). NLS1 is 
a subgroup of Seyfert 1 galaxies, but is characterized by large soft X-ray excess (see 
Fig.l in Boiler 2000). It has been proposed that black hole masses are systematically 
smaller, while luminosities are comparable to those of other Seyfert Is, leading to large 
L/L Edd in NLSls (Osterbrock & Pogge 1985; Boiler, Brandt, & Fink 1996; Mineshige 
et al. 2000). Our composite soft spectra may correspond for various spectra observed 
in NLSls. Furthermore, a small fraction of QSOs with optical-X-ray spectral index 
a Q x > 2 are observed (Yuan et al. 1998), which may be interpreted by our composite 
soft spectra. Comparison with GBHCs are another interesting, outstanding issue but 
is beyond the present investigation. 

Compared with the excellent work of Haardt & Maraschi (1993, hereafter HM93), 
our hard-state results are quite similar to theirs. This is easy to understand, since 
our study on the interaction between the disk and corona shows / ~ 1, which is a 
basic assumption in HM93. We also find soft-state spectra, where the corona is not 
coupled with the disk in the same way as that in HM93. Therefore, the spectra are 
different from HM93. Despite of the complexity, the basic differences can be briefly 
described as follows. HM93 consider the coupling between the disk and the corona 
through seed photons, which gives a relation between the soft radiation temperature 
Tp>, the corona temperature T, and the optical depth r, i.e. Fi(Tr,T,t) = C\. They 
find for 5eV < T R < 50eV and 0.01 < r < 1, the relation can be approximated as 
(160 2 + 40)r = 0.6 (where O = kT/m e c 2 ). For given T R and r, T is determined by this 
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approximation, and then the spectrum is obtained from the Monte Carlo simulations. 
By changing Tr and r within the limited range, they got somewhat similar spectral 
indices. In our model, we also consider the energy coupling between the disk and 
corona through the soft photons, i.e. Fi(Tr,T,t) = C\. In addition, we consider the 
coupling by the magnetic field and get the second relation, Fi(Tr, T, t) = C^M, M, f). 
Supplemented by the mass evaporation, we have the third relation, F 3 (T, r) = C3. 
Therefore, we uniquely determine T R , T, and r for given /, M, and M. Then, we 
perform the Monte Carlo simulation and obtain the spectrum and upward luminosity 
L(M, M, /). Using the fourth relation of that the luminosity equals to the accretion 
energy, i.e., L{M, M, f) = L G (M,M), we can also determine /. HM93 also take into 
account other factors like anisotropy of the seed photons. Here we concentrate on the 
underlying physics in order to determine the corona variables. 



5. Conclusion 

We improved our simple model presented in Paper I by introducing adjustable co- 
efficients for the energy density of the coupled seed photons and for the optical depth 
of plane-parallel corona. For given mass of black hole, accretion rate, and presumed 
coefficients, the model can determine the energy fraction dissipated in the corona and 
hence determine the structure of both the corona and disk. Then, we perform Monte 
Carlo simulations to trace the motion of photons in the corona and record the up- 
ward photons (emergent luminosity) and the backward photons. By comparing the 
presumed seed photon energy and the sum of intrinsic disk radiation and backward 
radiation energy, we adjust the coefficient of soft photon field and; by comparing the 
total accretion energy and the simulating emergent luminosity, we modify the effective 
optical depth in the slab geometry. Repeating this procedure, we finally obtain the 
self-consistent solutions for the corona and also for the disk, both of which are coupled 
by the magnetic field and the radiation field. 

We find two types of solutions for the disk and corona structure, corresponding 
for the hard-state and soft-state spectra. In the hard-state solution, the accretion disk 
is gas pressure-dominated. Most of the accretion energy is transferred to the corona 
by the magnetic field, and the intrinsic disk radiation is very weak. The corona and 
disk are tightly coupled through the radiation field and, consequently, the hard X- 
ray spectral indices are almost fixed at a ~ 1.1. In the soft-state solution, most of 
the accretion energy is dissipated in the disk and thus the disk is radiation pressure- 
dominated. The corona is weak, characterized by lower temperature and density, and is 
no longer strongly coupled with the disk through seed photons. The emergent spectral 
energy comes mainly from the un-scattered intrinsic disk radiation, peaking at UV and 
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soft X-rays. The soft states appear at only very high luminosity; while hard states 
appear at both high and low luminosities. A state between soft and hard, composed 
of an inner soft-sate solution and an outer hard-state solution, is also possible for a 
moderately luminous system. Our model predicts that the spectra is hard from low- 
luminosity objects with L < 0.2LEdd5 While the spectra can be either hard or soft 
from high-luminosity objects with L > 0.8LEdd- Spectra from objects with luminosity 
0.2L Edd < L < 0.8L Edd can be hard, and can also be moderately soft, which is softer 
at higher luminosity. 
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Fig. I. — Accretion-energy fractions dissipated in the corona, /, and in the disk, 1 — /, 
for M = 10 8 M Q and M = 0.1M Edd (corresponding to L = 0.07L Edd ). The disk is gas 
pressure-dominated. The figure shows that most of the accretion energy is dissipated 
in the corona. Only a very small fraction is dissipated in the disk. 
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Fig. 2. — Two solutions of / for M = 10 8 M Q and M = 2.0MEdd (corresponding to 
L = 1.4LEdd)- The lower panel shows that for a gas pressure-doimnated disk (hard 
state) a large fraction of accretion energy is dissipated in the corona; The upper panel 
shows that for a radiation pressure-doimnated disk (soft state) only a very small fraction 
of accretion energy is transported to the corona, while most of the energy is dissipated 
in the disk. 
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Fig. 3. — Coronal quantities along distance for M = 1O 8 M and M = O.lMEdd (corre- 
sponding to L — 0.07LEdd)- The corona is above a gas pressure-dominated disk with 
energy fraction channeled to the corona / ~ 1 shown in Figure 1. The coronal temper- 
ature is around 10 9 K and density around 10 9 cm~ 3 . The effective optical depth r* and 
Compton y-parameter y* are also shown, which hardly change with the luminosity. 
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Fig. 4. — Coronal quantities along distance for M = 1O 8 M and M = 2.0M E dd (cor- 
responding to L = 1.4L Edd ). The lower panel shows structure of the corona above a 
gas pressure-dominated disk with energy fraction channeled to the corona / ~ 1. The 
temperature and density are almost the same as that for small M shown in Figure 
3; while the Alfven speed is larger at higher M. The upper panel shows structure 
of the corona above a radiation pressure-dominated disk. Since most of the accretion 
energy is dissipated in the disk, the corona is weak with temperature ~ a few lO 8 ^ 
and density < 10 8 cm -3 in the inner region. In the outer region, T and n are relatively 
high due to a little larger fraction of local accretion energy is transfered to the corona. 



- 22 - 



46 I 1 1 1 1 1 1 i 1 i 1 1 1 1 1 r 




14 16 18 20 

Log(i//Hz) 



Fig. 5. — The hard-state spectra from a disk and corona around a black hole of 10 8 M Q 
at accretion rates M = 0.1M Edd (or L = 0.07L Edd ; Dotted curve); M = 0.5M Edd 
(or L = 0.35L E dd; Dashed curve); M = 1.0M E dd (or L = 0.7L Edd ; Solid curve); and 
M = 2.0M Ed d (or L = 1.4L Edd ; Dash-dotted curve). The hard X-ray emission is quite 
strong with a spectral index around 1.1. The spectral energy distributions are similar 
for different luminosities. 
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Fig. 6. — The soft-state spectrum (Solid curve) from a disk and corona around a black 
hole of 1O 8 M at accretion rate M = 2.0M E dd (or L = 1.4L E dd)- The spectral energy 
distributes dominantly in UV and soft X-ray since most of the accretion energy is 
dissipated in the disk. For comparison the corresponding hard-state spectrum is also 
shown in the figure (Dotted curve). Obviously, the X-ray spectra at soft state and hard 
state are of much difference. 
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Fig. 7. — The composed spectrum from a disk and corona around a black hole of 
10 8 M Q at accretion rate M = 0.5MEdd ( or L — 0.35LEdd)- The disk is radiation 
pressure-dominated in the inner region until 19i?s and gas pressure-dominated in the 
outer region, and hence the overlying corona is weak in the inner region and strong 
in the outer region. Dashed curve shows spectral energy distribution from the in- 
ner region, which is obviously contributed by the disk radiation; Dotted curve shows 
the contribution by the outer region, which is a typical radiation-coupled hard-state 
spectrum; The solid curve shows the total spectrum. The overall spectral energy distri- 
bution is between the hard sate and soft state, and the softness increases with accretion 
rate. 
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Fig. 8.— The hard-state spectrum for M = 1O 5 M and M = lM Edd (or L = 0.7L Edd ) 
(Dotted curve) compared to that for M = 1O 8 M at same accretion rate M = lM Edd 
(Solid curve). The spectral shapes are similar except that the disk radiation for M = 
10 5 M Q is at higher frequency due to higher disk temperature for lower black-hole mass. 
An example of £ = 20 Rs for M = 1O 8 M and M = lM E dd is also shown in the figure 
(Dashed curve), which overlaps in the standard case (£ = 10i?g). This indicates that the 
spectrum hardly changes with the thickness of corona since the hard-state temperature 
and optical depth given by our model depend on £ very weakly. 



